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' Abstract 

" ■ Quark-hadron phase transition in QCD in the presence of magnetic field is stud- 

, ied. It is shown that both the temperature of a phase transition and latent heat 

I decrease compared to the case of zero magnetic field. The phase diagram in the 

<~~^- plane temperature-magnetic field is presented. Critical point, T* = 104 MeV, 

Q^, = 600 MeV, for which the latent heat goes to zero, is found. 

' '\ PACS:11.10.Wx,11.15.Ha,12.38.Gc,12.38.Mh 

m ■ 

> : 

in \ 1 Introduction 

T— I ■ 

^ ■ QCD is essentially nonperturbative at low temperatures T < Tc {Tc is a quark-hadron 
2 ! phase transition temperature) and is characterized by confinement and spontaneous chiral 
00 ! symmetry breaking. In the hadron phase, at low temperature, the dominating contribu- 
^ I tion to the partition function of the system is given by lightest particles in the physical 
■ spectrum. In the case of QCD these particles are vr-mesons, which in the limit of massless 
^ ■ quarks are Goldstone excitations in the chiral condensate. It is a common method in the 
low temperature physics of QCD to use effective chiral theory [H El [3] , often called chiral 
perturbation theory (ChPT). 

One of the important questions is the phase structure of vacuum in the external 
magnetic field H. In the recent paper [4| it is argued that high magnetic fields eH ~ 
10^ -7- 10^ MeV^ are created in heavy ion collisions. Such magnetic fields should lead to 
effects that can be experimentally observed at RHIC. Also high magnetic fields eH ~ 
^QCD could exist in the early universe at the scale of strong interactions. Such high 
field strength can lead to new interesting phenomena at the stage of quark-hadron phase 
transition. At the same time it is interesting to study the influence of the external 
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magnetic field on the dynamics of strong interactions from purely theoretical point of 
view. Different nonperturbative phenomena in abelian magnetic fields were previously 
studied in[5l[6l[7l[8lSmM[ia[iara[l5l[l6l[I7l[T8l[iaEQ^ 

In this paper, we study the quark-hadron phase transition in QCD in abelian magnetic 
field. The physics of the considered phenomenon is the following. The plasma of hot 
quarks and gluons at T > Tc in the magnetic field is a thermodynamic system in a 
paramagnetic phase. On the other hand, at low temperature, T < T^, hadron matter, 
which mainly consists of scalar vr-mesons, is in diamagnetic phase. The paramagnetic 
phase is preferable thermodynamicly, because it minimizes the free energy (maximizes 
pressure). Therefore the temperature of the transition from hadron phase to the quark- 
gluon phase decreases as compared to the case of zero magnetic field H = 0. Thus, there 
is the analogy with the physics of condensed matter: confinement phase corresponds to 
the diamagnetic gas of scalar vr-mesons (we neglect the contribution of heavier hadrons), 
and deconfinement phase corresponds to the paramagnetic phase of quarks and gluons. 

2 Free energy of the QCD vacuum at T 7^ and H ^ 

The partition function of QCD in euclidian formulation in the presence of external abelian 
field can be written in a form (here T = is temperature) 

Z = exp 1^-^ dxi J d^xFl^ j[DB][Dq][Dq]exp\^- dxi j d^xC^ , (1) 
where QCD Lagrangian in the background field is 

^0 q=u,d 

here Qq is the charge matrix of quarks with flavor q = {u,d), and we omitted ghost 
and gauge fixing terms for simplicity. Free energy density is given by the expression 
(3VF{T,H,mg) = -InZ. 

Let us consider the hadron phase. At low temperature, T < T^, [Tc is the temper- 
ature of chiral phase transition) and at weak external field, H < nl^dr ~ the 
characteristic momentum in the vacuum loops is small, and the theory is described by 
the effective low-energy Lagrangian L^s [21 E] . 



Zeff[T,iJ] = e-'3^^^«[^'^] = Zo[iJ] [[DU]exp{- f dx, [ d'xL,f,[U, A]}. 

J Jo Jv 



(3) 



Lefj can be represented as a decomposition in a series of powers of momentum (deriva- 
tives) and masses 

Leff = L(2) + L(^)+L(6) + ... (4) 
The leading term in (jl]) is a Lagrangian of nonlinear a-model in external field 

1(2) ^ ^Tr(V^t/+V^f/) + Si?eTr(Mf/+), 
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V^U = d^U-z[U,V^]. 



(5) 



Here U is unitary SU (2) matrix, F.,^ = 93 MeV is a pion decay constant, and the parameter 
S has the meaning of quark condensate, S = | (uu) \ = \ (dd) \ . External abehan magnetic 
field H, directed along z axis, corresponds to V^(x) = {t^ /2)eA^{x), and vector-potential 

is taken in the form A^{x) = 5^2Hxi. We will furhter neglect the breaking of isospin 
symmetry of strong interactions and consider masses of light u- and d- quarks equal, 
TTiu = TTtd = rriq, thus the mass matrix diagonal, M = iriql. 

In one-loop approximation it is sufficient to use the decomposition of L^a up to 
quadratic over pion field terms. In exponential parametrization of the matrix U{x) = 
exp{iT°"n'°'{x) / Ftt} we find that 



(6) 



where we introduced the fields of charged vr^ and neutral 7r° mesons 



(7) 



The QCD partition function ([T]) in one-loop approximation of the effective chiral theory 
takes the fornu. 

Z^^[T,H] = Zp^Zo[H] [ [Dn°][Dn+][Dn-]exp{- f dx^ [ d^xL^^) [vr, A]}. (8) 

J Jo Jv 

Here the partition function is normalized to the case of perturbation theory at T = 0, 
H = 

Zpr=[deti-dl + M'^)]-'/\ (9) 
Integrating ([8]) over vr- fields one gets 



Z^s = Zp'^Zo[H][detT{-dl + M'^)]-'/'[detT{-W + ^')]-^ 



(10) 



where = — ieA^ is a covariant derivative, and the subscript T means that the 
determinant is evaluated at finite temperature T using standard Matsubara rules. Using 
and regrouping multipliers in fllUp we arrive at the following expression for Z^ 

-1/2 r 1 w I 19 , ,^9Nn-l 



Z^^[T,H] = Zo[H] 



detri-dl + Ml) 
dei{-dl + Ml) 



det{-\Dy + Ml] 
dei{-dl + Ml) J 



X 



detT( — 




? + Ml) 


det(-| 




2 + Ml) 



Free energy then takes the form [7] 



F^^iT,H) 



+ F^o{T) + F^4H) + FUT,H). 



(12) 



^The partition function describes charged tt^ and neutral 7r° Bose gases in magnetic field 
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Here F^o is the free energy of massive scalar boson 



FAT)=T J ^ln(l-exp(-vWM|/T)), 



(13) 



is Schwinger's result for the vacuum energy density of charged scalar particles in 
magnetic field 

eHs 



- 1 



(14) 



sinh(eifs) 

Next, it is technically not hard to generalize the case of the vacuum H = 0, T = for 
charged vr^-mesons 

F = Trln(p2 + ^2(p)) (15) 

to the case of H ^ 0,T ^ 0. Omitting the details, we will note that it corresponds to the 
following substitutions 

P4 — >■ cufc = 27ikT, {k = 0, ±1, ...), 
LUo = Vp^Tm^ ^„ = ^M+JI^T^Hi2^^TT) (16) 

and 

^ 27r ^ ^ J ^ 27r ' 

n=0 k=-oo 

where the degeneracy multiplicity of Landau levels eH/2iT is taken into account. Summing 
over Matsubara frequencies, one obtains the following result for the diamagnetic part of 
free energy of charged Bose gas 

TTrp ^ /"OO 

FUT,H) = ^ V / dA;ln(l -exp(-a;„/T)), 

^ n=0 



uJn = VP + M2 + ei7(2n + l), (17) 

here Un are Landau levels of vr^-mesons in a constant magnetic field H. 

Expanding ln(...) in the integrand of ( fT3l) . ( |T7I) in a series, one gets the following 
expressions: 



n=l 

and 



^dia = - V 5^ + eH{2n + 1) J] ^^^i ^ + eiJ(2n + 1) 

n=0 fc=l ^ ^ 



(19) 



where K„ is the Macdonald function. 
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3 Quark-hadron phase transition in magnetic field 



In order to determine the temperature of phase transition we will write down the pressure 
in two phases. At zero chemical potential pressure is equal to minus free energy. 

In the confinement phase pressure of vr-mesons in magnetic field takes the form (see 

m) 

P^{T, H) = PAT) + P.±{H) + PdiaiT, H), (20) 

where the pressure of neutral gas of 7r°- mesons is P.,^o{T) = —Ft^o[T). Renormalized 
contribution of vacuum polarization of charged vr^-mesons to the pressure, which does 
not depend on temperature (Schwinger polarization) is 



IGtt^ In s 



1 + -{eHsf 



sinh (eifs) 6 



(21) 



Diamagnetic term in the pressure, which comes from charged tt^ mesons is P^iaiT, H) = 

— Fdia{T, H). 

Fermion (quark) determinant at finite temperature in the magnetic field can be con- 
sidered in a similar way. Then the pressure in the quark-gluon plasma phase has the 
form 

P,i{T, H) = 2{Nl - 1)^T^ + P,{H) + P,,,„(T, H). (22) 

The term proportional to oc T'^ comes from the gas of hot gluons, Pq{H) is the contribution 
from vacuum polarization of quarks to the pressure 



q=u,d 

and Ppara is a paramagnetic term 



°° ds 



S3 



\eg\Hscth{\eq\Hs)-l-^{eqHsY 



(23) 



p..^2iv.5:^Tf:5:/;3in(i.e- 



^^o.tTi'^- ' (24) 



ujg = pl + ml + \eq\H{2n + 1 + 2a). 

An important physical property of the phase transition is the rearrangement of the 
nonperturbative QCD vacuum. Due to the scale anomaly in the trace of the energy- 
momentum tensor the new dimensional quantity, gluon condensate {G^) = ((s'G^y)^), 
emerges in QCD. Nonperturbative energy density of the vacuum is related to (G^): 

^ {G'), (25) 



1287r 



where b = {3Nc — 2Nf)/3, iVc = 3 is the number of colors, and b = 29/3 for Nf = 2. 

Energy density of vacuum is the negative quantity, and the state with (G^) 7^ turns 
out to be thermodynamically advantageous. Theoretical studies [211 122] and numerical 
computations in the lattice approximation of QCD [23] show that at the point of phase 
transition Tc one part of the condensate (chromoelectric part) turns to zero, while the 
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chromomagnetic condensate remains almost unchanged as compared to the case T = 0. 
In the vacuum at T = {{E°;Y) = {{H"')'^), and therefore vacuum energy density above 
phase transition appears to be less in magnitude than below phase transition, and the 
difference is approximately 

Taking this into account, the quantity —Asy should be added to the equation of state 
in the plasma phase. Thus, the pressure in the quark-gluon plasma state is given by the 
expression 

P2{T) = Ppi{T,H)-Aey. (27) 

Equation ( 1271) is similar to the equation for the phase transition in MIT bag model, where 
bag constant B plays the role of Asy. 

Phase transition temperature can be found from the condition of equality of pressures 
in both phases 

P^{T,,H)=P2{T,,H). (28) 

Let us now consider "weak" magnetic field, eH -C T^?. Then the Schwinger contri- 
bution to the pressure can be neglected. In the weak field one can use Euler-Maclaurin 
formula for P^ia 



-F{a) + ^ F(a + n) ^ / dxF{x) - -F\a 

^ n=l 

and (|T71) can be rewritten in the form 

Jo Vx^ + z^{e^^^ -1)' 
For the paramagnetic term Ppara we find in the weak field 

P..„-Po.(T,if^O) + ^^^^^/,(^ 



) (29) 



(30) 



dx 



(31) 



Va;2 + z2(ev^^^ + l) 



where 1^^ = = m and = (e^ + e^)/e^ = (| + |) = |, and pressure Poq(T, H = 0) is 
given by □ 

Po.(T, H = 0) = x^dx In (l + e'-^/-) , ^^^^ 



^Using x^dx\^{l + e-'') = g^, one finds in the limit m = that Poq{T,H ^ 0) = ANc—T"^. 
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In the absence of magnetic field and in the chiral hmit one finds pi] 



(7-3)(7r790) 



1/4 



(33) 



Here 7 = 2- (A*"^ — 1) + (7/8) - 2 ■ 2-Nc-Nf is the number of independent degrees of freedom 
of quarks and gluons, and 7 = 37 for Nc = 3, Nf = 2. Lattice calculations give the value 
(G^) = 0.87 GeV^, and one finds from (133!) phase transition temperature ~ 177 MeV 
Hi H = 0. 

The influence of magnetic field can be taken into account in the first approximation 
by redefining vacuum energy density 



where 



V 



As^ = As, - {eHyv, 



127r2 



m 



(34) 



(35) 



V = Q.l- 10^2 for = UOMeV and m ^ 5MeV. 

Considering the term {eH^V as a perturbation, one finds from (155]) and (131|) that 
the relative shift of the deconfinement phase transition temperature in the presence of 
magnetic field is 



T, 



H 



Tr 1- 



4Ae, 



(36) 



and V/AAe^, - 9.2GeV-\ 

Thus, we see that the presence of a magnetic field leads to a decrease of the quark- 
hadron phase transition temperature, and AT/Tc ~ 10~^(eif)^/Ae^. 



4 Results of numerical simulations 

Equations (120!) and (1271) allow to evaluate the pressure in both phases, and to numerically 
find the dependence of phase transition temperature on the magnetic field from (1281) . 
Results of numerical calculations for = 3, Nf = 2 are presented in Fig. HJ Phase 
transition temperature, as discussed above, decreases with increasing external magnetic 
field. 

Thermodynamics in each phase is defined by the pressure, and we can evaluate energy 
density and latent heat in both phases. Energy density is given by 



dP2 

62 = T - P2 

dT 



(37) 



Latent heat equals to the difference of energy densities of two phases in the point of phase 
transition: 

Ae(i7)= (62-^1) It.(/.) (38) 
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Figure 1: Quark-hadron phase transition temperature vs external magnetic field. Critical 
point, T^, = 104 MeV, \feHl = 600 MeV, where latent heat turns to zero, is marked 
with the dot. Solid line corresponds to the first order phase transition, and dashed line 
corresponds to the crossover. 

The dependence of Ae{H) is plotted in Fig. [2J The value of the magnetic field 
y/eH^ = 600 MeV where latent heat turns to zero corresponds to the critical point, at 
which first order phase transition changes to the crossover. 

5 Conclusion 

We have studied the quark-hadron phase transition in QCD in the presence of external 
magnetic field, and have shown, that the temperature of the phase transition decreases in 
comparison to the case of zero magnetic field. Equation fl28l) was solved numerically, the 
phase diagram in the plane temperature-magnetic field and critical point were found. 

As was shown above, there are two phases in the presence of magnetic field: diamag- 
netic phase below Tc and paramagnetic above Tc. Correspondingly magnetic susceptibility, 
X = —d^P/dH"^, changes it's sign at the critical temperature. Thus, magnetic suscep- 
tibility may be considered as the order parameter of the model of thermal QCD in the 
presence of magnetic field. 

It is known from lattice calculations that there is a crossover for finite temperature 
QCD with physical quark masses. In the presence of magnetic field there are additional 
magnetic terms in the pressure, which give different contribution to the energy density in 
two phases. Thus we expect that a crossover is replaced by a first order phase transition. 
Analogous phenomenon was found in [24j, where it was shown that chiral transition 
changes from crossover to the weak first order transition in the linear sigma model in a 
magnetic field. 
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Figure 2: Latent heat Ae{H) vs external magnetic field V eH 

The following remark should be made. From |7j it is known that the chiral phase 
transition temperature grows with the magnetic field. As it follows from lattice calcula- 
tions, deconfinement and chiral phase transitions take place at the same temperature Tc 
in case of zero magnetic field H = 0. On the other hand, as was shown above, quark- 
hadron phase transition temperature at nonzero magnetic field is lower than in case of 
H = 0. Thus, chiral and quark-hadron phase transition temperatures are separated in 
the presence of magnetic field. Correspondingly, there appears a temperature interval, 
where the quark-hadron phase transition is already passed, but the chiral symmetry is 
still broken. This phenomenon may be important for the consideration of quark-hadron 
phase transition in heavy ion collisions and in the early universe. 
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